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In an electron system coupled with anharmonic phonons, i.e., rattling, inverse isotope effect on the 
Kondo temperature Tk is found to occur by the numerical evaluation of the Sommerfeld constant 7 of the 
Anderson-Holstein model. For the anharmonic potential of an oscillator with mass AI in which large 7 
has been found to be almost independent of an applied magnetic field, 7 is significantly suppressed when 
M is increased, i.e., Tk is enhanced due to the relation of J'~^Tj^^ in the Kondo problem, leading to the 
inverse isotope effect on Tk. Since this phenomenon does not occur for harmonic phonons, it can be a key 
experiment to prove the relevance of rattling to magnetically robust heavy electron state. 
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In traditional heavy-fermion materials such as CeCug,' 
the electronic specific heat coefficient 7 is significantly sup- 
pressed by an applied magnetic field. This can be understood 
within the standard Kondo effect,^ i.e., the singlet formation 
of local magnetic moment due to the coupling with conduc- 
tion electrons.^ From the viewpoint of the standard Kondo 
effect originating from spin degree of freedom, it is quite un- 
usual that heavy effective mass is independent of the applied 
magnetic field. Thus, the magnetically robust heavy-fermion 
phenomenon observed in SmOs4Sbi2'^ has attracted much at- 
tention of condensed-matter theorists due to the renewed in- 
terest on the realization of non-magnetic Kondo effect. 

The theoretical research of the Kondo effect has a history 
longer than forty years. ^ It has been widely recognized that 
the Kondo-like phenomenon generally occurs in a conduction 
electron system hybridized with a localized entity with inter- 
nal degrees of freedom. For instance, when electrons are cou- 
pled with Einstein phonons, a double-well potential is formed 
in an adiabatic approximation, leading naturally to a two-level 
system, in which Kondo has first considered a possibility of 
non-magnetic Kondo behavior''-^ This two-level Kondo sys- 
tem has been shown to exhibit the same behavior as the mag- 
netic Kondo effect. '^•'^ In the context of heavy-fermion behav- 
ior in A15 compounds, Kondo effect with phonon origin has 
been discussed.'" " Further efforts to include the low-lying 
levels of local phonon has been done to understand the quasi- 
Kondo behavior in electron-phonon systems.'^ 

Due to the discovery of the magnetically robust heavy- 
fermion phenomenon in SmOs4Sbi2, the research on the 
Kondo effect with phonon origin seems to revive in recent 
years. This phenomenon occurs on the filled skutterudite 
structure, in which rare-earth guest atom is surrounded by the 
cage composed of twelve pnictogens. Then, the rare-earth ion 
easily moves around potential minima or wide bottom of the 
potential inside the pnictogen cage. This anharmonic motion 
is called rattling, which is believed to have significant effect 
on electronic properties of filled skutterudites. 

As an extension of the two-level Kondo problem, the four- 
and six-level Kondo systems have been analyzed to clar- 
ify the effect of rattling in filled skutterudite.'^ ''* The ori- 
gin of heavy mass has been discussed on the basis of the 
Anderson-Holstein model. '^ '"^ The present author has dis- 
cussed the magnetic field dependence of Sommerfeld constant 



in the Anderson-Holstein model. Then, it has been shown 
that large 7 becomes magnetically robust when the bottom of 
the potential becomes flat due to the effect of anharmonicity. 
From these efforts, it has been gradually understood that mag- 
netically robust heavy-fermion phenomenon actually occurs 
in electron-rattling systems, but we do not have an evidence 
for the relevance of rattling to this intriguing phenomenon. 
Thus, it is meaningful to consider an experiment which con- 
firms the contribution of rattling to the Kondo effect. 

In this letter, we show that inverse isotope effect occurs in 
the magnetically robust heavy electron state. For the purpose, 
we numerically evaluate the Sommerfeld constant 7 of the 
Anderson model coupled with local oscillation of guest atom. 
For the anharmonic potential in which 7 is almost indepen- 
dent of an applied magnetic field, when we increase the mass 
of the oscillator M, we observe the decrease of 7, i.e., the 
increase of the Kondo temperature Tk due to the relation of 
7~Tj^^ in the Kondo problem. Namely, in the expression of 
TkocM^^, the exponent rj becomes negative. This is the in- 
verse isotope effect on Tk, which can be an evidence for the 
relevance of anharmonic phonons to the Kondo effect, since 
it does not occur for harmonic phonons. 

Let us first discuss how the isotope effect appears in the 
local electron-rattling state. The local term is expressed as 
Hioc=Un^ni + gnQ + P'^/{2M) + V{Q), where U denotes 
Coulomb interaction, ncr=d'ld„, d„ is an annihilation oper- 
ator of localized electron with spin a at an impurity site, 
g is the electron-phonon coupling constant, n=n^-\-ni, Q is 
normal coordinate of the oscillator, P is the coiTesponding 
canonical momentum, M is the mass of the oscillator, and 
V{Q) is the potential for the oscillator Here we express V{Q) 
as V{Q)=kQ^ /2 + k^Q'^ + kgQ^, where A; is a spring con- 
stant and we consider fourth- and sixth-order anharmonicity, 
fc4 and fcg, in the potential. Note that we use such units as 
h=kB=l in this paper 

By following the standard procedure of quantization of 
phonons, we introduce the phonon operator a defined through 
Q={a + a^)/y/2ujM, where uj is the phonon energy, given by 
UJ= y/kjM. Then, we obtain 

^^loc = Unini + uj\/a(a + a^)n + Lu{a^a + 1/2) 

(1) 

where a indicates the non-dimensional electron-phonon 
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coupling constant defined by 0=0^ j (2Mlj^') and non- 
dimensional anharmonicity parameters, (i^ and /Sg, are given 
by ^4=^:4/ (4M2cj3) and liQ=kfi j (%M^ uj% respectively. 

Here we explain M dependence of parameters. First we 
note the relation lo cx M^^/^, since we usually consider that 
the spring constant fc does not depend on M. From the defini- 
tion, a itself depends on M as a oc 

Ml/2, 

auj does not 

depend on M. Here it is instructive to recall the isotope effect 
in BCS superconductor Namely, the M dependence of su- 
perconducting transition temperature occurs only through 
Lo, since the effective attraction between electrons mediated 
by harmonic phonons is given by auj. Then, we obtain the 
famous formula oc M^^/^ for the BCS superconduc- 
tor Concerning anharmonicity parameters, we also consider 
that fc4 and fcg are independent of M. Thus, we obtain that 
Pi cx and /Jg oc Af-^. 

In this paper, we set a)=cj/-^/m, a=a-Jrn, [3/i=(3i/ \pm, and 
PQ=l36/m, where m denotes the mass ratio of the oscillator 
when we substitute the oscillator atom with its isotope. The 
quantities without tilde denote those for m=l. For the actual 
calculations, we define the phonon basis as | £) = (a^^ ) ^ 1 0) / \/T\, 
where I is the phonon number and |0) is the vacuum state. 
The phonon basis is truncated at a finite number, which is set 
as 1000 in this paper. 

We briefly explain the change of the potential shape due to 
j3i}^ The potential is rewritten as V{q)= aijj{q^ + WaP^q'^ + 
64:0^ Peq^), where q is the non-dimensional length given 
by q=Quj^/g. Note that V{q) is independent of M, since 
auj=auj, a(3i=a(3i, and a"^ l3Q=a'^ (3^. In Fig. 1(a), we show 
V{q)/{auj) for several values of Pi for a=2 and (3q=1Q^^. For 
(3^=0, V{q) has a single minimum at q=Q. When we decrease 
Pi, shoulder-like structure begins to appear around q=±4 
and the bottom of the potential becomes relatively wide. For 
/34<— 0.00274, potential minima at q^O appear When Pi is 
smaller than —0.003, a couple of minima at q^Q are grad- 
ually deep. In the following, the potential in the region of 
-0.003</34<-0.002 is called the rattling type. 

As mentioned above, there occurs attractive interaction C/ph 
between electrons mediated by phonons. The effective inter- 



action Ucs is evaluated by Ucff=U~ 

,(0) 



-[/ph=£'o*'^- 



-E. 



(0) 



(0) 



where En is the ground-state energy of iJioc for local elec- 
tron number n. For the case of harmonic potential (Pi=Pe=0), 
C/ph is analytically evaluated as t/ph=2acj. Thus, we obtain 
UcS=U — 2aw, which does not depend on M, when we sim- 
ply assume that U is not affected by M. 

On the other hand, for the anharmonic potential, significant 
M dependence appears in Ucs, which will be a source of the 
isotope effect on the Kondo temperature. In order to visualize 
the M dependence of Ucs, we consider the variation of C/eff 
due to the increase of M, given by (5f/off=C4ff — t/off (1), 
where Ucfiim) is the effective interaction for m. Since U is 
assumed to be unchanged by M, the variation occurs only 
through the change of [/ph. 

In Fig. 1(b), we plot SUcS as a function of Pi for m=1.05, 
Pq=1Q^^, a=2, and U /uj=\Q. For harmonic potential, 5Ucs is 
zero, as shown by the red line. For anharmonic potential, Ucs 
is decreased, i.e., [/ph is increased, due to the increase of M 
for the potential with single minimum, while L/off is increased 
with the increase of M for the potential with a couple of deep 
minima at q 7^ 0. For the rattling-type potential, the variation 
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Fig. 1. (Color online) (a) Anharmonic potential vs. q for .several values of 
f3i. Other parameters are /36=10~^ and a=2. (b) Variation of the effective 
interaction SU^ft as a function of f3i for m=1.05, /3g=10~^, a=2, and 
U /w=lO. The red line denotes the result for harmonic potential. 



of C/ph due to the increase of M is relatively large, since the 
flat potential easily deforms by electron-phonon coupling. 

Let us now consider the hybridization between localized 
and conduction electrons. The Hamiltonian is called the 
Anderson-Holstein model, given by 

J 



-fin- 



loC; 



(2) 



where Sk denotes the dispersion of conduction electron, Cka 
is an annihilation operator of conduction electron with mo- 
mentum k and spin cr, /i is a chemical potential, and V is the 
hybridization between conduction and localized electrons. In 
the following, we set V / D=Q.25, where Z? is a half of the 
conduction bandwidth. Note that we adjust to consider the 
half-filling case, even though we do not mention explicitly. 

For the evaluation of the Sommerfeld constant 7 of the 
Anderson-Holstein model, we employ a numerical renormal- 
ization group (NRG) method,^' where we include efficiently 
the conduction electron states near the Fermi energy by dis- 
cretizing momentum space logarithmically. In the actual cal- 
culation, we introduce a cut-off A for the logarithmic dis- 
cretization of the conduction band. Due to the limitation 
of computer resources, we usually keep only L low-energy 
states. In this paper, we set A=5 and L=2000. By using the 
NRG method, we calculate the Sommerfeld constant 7 at a 
low temperature T. Note that T is defined as rjv=A"'^"i^/2 
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Fig. 2. (Color online) (a) Ratio of the Sommerfeld constant ■y{Tn) /■y(l) 
vs. /34 for U/D=2, a=2, uj/D=0.2, /36=10-^, and m=1.05. Numerical 
results are shown by red circles. Solid curve is obtained from the analytic 
expression of Tk of the effective s-d model, except for the region with 
UeB < 0. (b) Exponent rj in the relation of Tk oc M"'' vs. /34 for several 
values of U. Other parameters are a=2, aj/D=0.2, and /36=10~^. 



in the NRG calculation, where N denotes the number of the 
renormalization step. We actually evaluate 7 through the rela- 
tion of j=C/T at a low temperature T, where C is the specific 
heat of localized electron. Namely, it is necessary to obtain C 
with high accuracy. For the purpose, we evaluate it by the 
numerical derivative of the entropy S through the relation of 
C=dS/d\ogT={SN-i - S'N)/log(Ai/2), where Sn is the 
entropy at the step N. 

Here we summarize the previous NRG results of 7 for 
771=1.^" We have found that phonon-mediated attraction is 
largely enhanced for -0.0027</34<-0.0021, in which the 
potential has wide bottom and SUcff is negative with large 
absolute value. There occurs a strong cancellation between 
Coulomb repulsion and the phonon-mediated attraction, just 
when the potential shape is deformed from the rattling type. 
We have shown that in such a situation, spin and charge fluctu- 
ations are comparable to each other, leading to the emergence 
of electron-rattling complex state with large and magnetically 
robust 7. In the present notation, magnetically robust large 7 
has been found ai'ound at /34«-0.0021 and -0.0027. 

Now we consider the isotope effect on 7 for rat- 
tUng phonons. In Fig. 2(a), we show (34 dependence 
of 7(m)/7(l) for U/D=2, a=2, uj/D=0.2, and m=1.05, 
where 7(m) is the Sommerfeld constant for m. For 



/34>-0.0021 and /34<-0.0027, we observe 7(to)/7(1)>1. 
For -0.0021 >/34>-0. 0027 with relatively wide bottom in 
the potential, when m is increased, 7 is decreased, i.e., Tk is 
increased due to the relation of "f^T-^^ in the Kondo prob- 
lem. This is the inverse isotope effect, since in the standard 
isotope effect, the relevant temperature is decreased with the 
increase of M, as observed in the BCS superconductor. 

The inverse isotope effect is caused by the further de- 
crease of negative effective interaction Ucb in the region of 
-0.0021>/34>-0. 0027 due to the increase of M, as found in 
Fig. 1(b). Note that /34>-0.0021, SUcS is negative, but UeS 
is positive with large absolute value. In this case, the change 
in Tk is determined by the exchange interaction, not by UcS- 
When Ucff is positive and is larger than the width of the vir- 
tual bound state, the Anderson-Holstein model is reduced to 
the isotropic s-d model with the exchange interaction J,'^ ex- 
pressed as 



E. 



(0) 



(3) 



where \^n^) is the £-th eigenstate of iJioc for electron number 
n with the eigenenergy En ■ Then, Tk is given by the well- 
known formula Tk = De"^/^^''"'-', where p is the density of 
states at the Fermi energy. The solid curve in Fig. 2(a) indi- 
cates [Tk(to) /Tk(1)]^^. Except for the region with negative 
[/off , the numerical results are reproduced by the rough esti- 
mation of Tk- For C/eff<0, analytic expression of Tk is not 
obtained, but we observe the tendency in the solid curve to- 
ward the decrease of [Tk(to)/Tk(1)]^^. Thus, it is believed 
that the inverse isotope effect is the essential feature of the 
anharmonic potential. 

In Fig. 2(b), we show the dependence of the expo- 
nent 77 in the formula of the isotope effect Tk oc M^^. 
For the evaluation of rj, it is convenient to use the rela- 
tion 77=log[7(m) /7(l)]/ logTO for m near unity. Here we set 
m=1.05. When 7(m)/7(l) is smaller than unity, the exponent 
becomes negative, directly suggesting the inverse isotope ef- 
fect. When we increase U, the region with negative 77 is not so 
rapidly suppressed, as found in the results for U/D=2.2 and 
2.5. In particular, for U/ D=3, Ucb is found to be always posi- 
tive in the present parameters, but we still find the negative rj. 
Thus, the effective attraction is not the only condition for the 
appearance of the inverse isotope effect. Rather, the key issue 
is the heavy electron state dressed by rattling phonons, when 
the potential bottom becomes wide. From the above results, 
77 is negative only for the rattling-type potential, including the 
region with magnetically robust 7. 

In order to confirm that the inverse isotope effect occurs 
only for anharmonic phonons, we discuss the results for har- 
monic phonons, i.e., /34=/36 =0. Roughly speaking, the elec- 
tronic state is understood by the comparison of U and 2auj. 
For U > 2auj, the local electron state has double degeneracy 
of spin degree of freedom, even though the Coulomb interac- 
tion is weakened by the attraction mediated by phonons. In 
this region, the model is described by the s-d model with the 
exchange interaction J. On the other hand, for U < 2auj, 
the vacant and double-occupied states are degenerate, lead- 
ing to the situation of the two-level Kondo system. In such 
a case, the Anderson-Holstein model is effectively reduced to 
the anisotropic s-d model, in which the longitudinal exchange 
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Fig. 3. 7(m)/7(l) vs. U for /34=/36 =0, m=1.05, a=2, and a;/Z)=0.2. 
Solid curve is estimated from Tk of the effective s-d model." 



interaction Ji is different from the transverse component J2. 
Due to the immobihty of bi-polaron in comparison with sin- 
gle polaron, J2 is smaller than Ji. As for more details, readers 
can consult with Ref. 19. 

In Fig. 3, we depict the U dependence of '-f{m)/j{l) for 
m=1.05, a=2, uj/D=0.2, and /34=Pe=0. For comparison, we 
depict the curve obtained from the analytic forms of Tk of 
the effective s-d models.'^ Since the expressions of Tk are 
obtained in the strong-coupling limit, the numerical results 
should not perfectly agree with the analytic ones. Neverthe- 
less, the overall behavior is well reproduced by the solid 
curve, even though the present situation is not in the strong- 
coupling limit. It is not surprising that the isotope effect is also 
found in Tk for harmonic phonons, since the virtual phonon 
excitations are included in the expressions for J, Ji, and 
J2. Here we emphasize that '~f{m)/j{l) is always increased, 
namely, Tk(to)/Tk(1) is decreased, when we increase m. 
Thus, the isotope effect on Tk appears for harmonic phonons, 
but Tk is always decreased with the increase of m. 

In order to detect the inverse isotope effect, we should es- 
timate experimentally 7 from the specific heat measurement, 
e.g., of SmOs4Sbi2 with isotope of Sm, since we consider the 
oscillation of guest atom. The idea is simple, but it is challeng- 
ing to find the direct evidence of the relevance of rattling to 
the Kondo effect. In the present paper, we have considered the 
isotope effect on 7, but we have also checked that the same ef- 
fect appears in magnetic susceptibility x, although the results 
are not shown here. Since the correspondence between x and 
Tk is clear in general, it may be more practical to discuss ex- 
perimentally the isotope effect on x- In any case, we believe 
that the change of 7 and x in the compound including isotope 
is detectable in experiments. 

Three comments are in order (i) We should pay due at- 
tention to the effect of local interaction, when we discuss the 
actual periodic system from the impurity model. For instance, 
the magnitude of 77 in the periodic system may not be so large 
in comparison with the present results. The point is that rj be- 
comes negative for magnetically robust heavy electron state, 
(ii) When we substitute Sm with its isotope in SmOs4Sbi2, 
Tk may be distributed, which makes it difficult to detect the 



isotope effect on 7 and x- This point seems to be related to 
the improvement of sample quality due to the increase of fill- 
ing fraction of rare earth atom in the pnictogen cage of filled 
skutterudites.^^ (iii) In general, Tk should be distinguished 
from the so-called coherence temperature concerning the for- 
mation of heavy electron state. The relation between Tk and 
observables is not so clear in actual heavy-electron materials. 
When rattling is relevant to the heavy electron state, the iso- 
tope effect may significantly appear in the coherence temper- 
ature. This point may be discussed on the basis of the periodic 
Anderson-Holstein model, but it is one of future issues. 

In summary, the Sommerfeld constant 7 of the Anderson- 
Holstein model has been evaluated by the NRG method. For 
the rattling-type potential of the oscillator, when we increase 
the mass of the oscillator, we have found the decrease of 7, 
i.e., the increase of Tk. The same effect also occurs in the 
magnetic susceptibility x- Then, we have concluded that the 
inverse isotope effect on Tk occurs, if the magnetically robust 
large 7 originates from rattling phonons. We expect that the 
measurements of 7 and x in SmOs4Sbi2 including Sm isotope 
will be performed in future. 
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